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Abstract 



fvl I In this paper, we propose a globally hyperbolic regularization to the general Grad's 

moment system in multi-dimensional spaces. Systems with moments up to an arbitrary 
order are studied. The characteristic speeds of the regularized moment system can be 
analytically given and only depend on the macroscopic velocity and the temperature. 
The structure of the eigenvalues and eigenvectors of the coefficient matrix is fully 
clarified. The regularization together with the properties of the resulting moment 
C^ ' systems is consistent with the simple one-dimensional case discussed in [1]. Besides, 

all characteristic waves are proven to be genuinely nonlinear or linearly degenerate, 
and the studies on the properties of rarefaction waves, contact discontinuities and 
shock waves are included. 

> 

\j^ . Keywords: Grad's moment system; regularization; global hyperbolicity; characteris- 

es^ ' tic wave 

m 
p 
en ' 1 Introduction 

^^ ' The kinetic gas theory, which is based on the Boltzmann equation, is one of the funda- 

mental tools in modelling non-equilibrium processes. Nevertheless, in most cases, a direct 
numerical discretization of the Boltzmann equation leads to unacceptable computational 

^ , costs. In 1940s, Grad [7] proposed the moment approximation of the distribution func- 

tion, trying to establish a series of intermediate models between the fluid dynamics and 
the kinetic theory. However, due to a number of defects in Grad's 13-moment equations, 
such as the appearance of unphysical subshocks, nonexistence of an entropy function, and 
lack of global hyperbolicity, not much attention is paid to the moment method in the last 
century. 

In the recent twenty years, as the investigation into the moment method becomes 
deeper, various "regularizations" are proposed to challenge the traditional accusations on 
the moment method. A list of relevant publications can be found in the references of 
|13j . Recently, we are interested in the large moment system together with its numerical 
methods [21 [5l |H [3], and it is found that the lack of the well-posedness due to the loss 
of global hyperbolicity is a major obstacle in our simulations, especially for large Mach 
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number gas flows [5]. Torrilhon [T3] provided a 13-moment hyperbolic moment system 
based on multi-variate Pearson- IV distributions, but it seems unlikely to extend the same 
technique to systems with large number of moments. As discussed in [13], Levermore [9] 
gave a partial answer to the question of hyperbolicity of large moment system based on 
a maximal entropy distribution function. However, the analytical forms of Levermore's 
equations cannot be obtained once the number of moments is greater than 10. While 
exploring the method ensuring the hyperbolicity of the moment system, we discovered [T] 
that the structure of the characteristic polynomial of Grad's moment equations with one- 
dimensional microscopic velocity is rather simple; thus a globally hyperbolic regularization 
can be achieved by simply adding two terms to the equation of the highest order moment. 

In this paper, the results in [1] are extended to the multi-dimensional space. For 
multi-dimensional moment systems, the regularization method is consistent with the one- 
dimensional case. Due to the complexity of the moment systems, this paper is mainly 
devoted to a rigorous proof of the hyperbolicity of regularized moment system for any 
space dimensions and an arbitrary order of moments. The result is obtained by firstly 
restricting the spatial variable in the one-dimensional space, and then it is generalized 
to the multi-dimensional space using the rotation invariance of the regularized system. 
For the case of one-dimensional spatial variable, the structure of the coefficient matrix 
is similar as the Hessenberg matrix, which enables us to calculate the eigenvectors for a 
given eigenvalue. Then, the hyperbolicity of the moment system follows by counting the 
number of linearly independent eigenvectors. At the same time, the expressions of all 
characteristic speeds are obtained, each of which is a sum of the macroscopic velocity and 
the square root of the temperature scaled by a zero of the Hermite polynomial. Besides, we 
prove that each characteristic filed of the hyperbolic moment systems is either genuinely 
nonlinear or linearly degenerate, and some properties of the rarefaction waves, the contact 
discontinuities and the shock waves are investigated. 

The rest of this paper is arranged as follows: in Section [21 a brief review on the 
moment methods of Boltzmann equation and the results in [1] are presented. Section 
[3| gives the globally hyperbolic regularization for moment system with one-dimensional 
spatial variable and multi-dimensional microscopic velocities. And in Section [U the result 
for full multi-dimensional moment system is proved. The study on the characteristic waves 
is carried out in Section [5) 

2 Preliminaries 

In this section, a concise introduction of the Boltzmann equation is presented. And then 
some results of the work on the moment method in [5l [1] are briefly reviewed. 

2.1 Moment methods for Boltzmann equation 

Let the motion of particles be depicted by the distribution function /(t, x, ^) governed by 
the Boltzmann transport equation 

^ + E^i^ = ^(/'/)' *^^^' ^^^^^""^ (2-1) 

where t denotes the time, x = (xi,--- ,xd) and $, = (^i,--- ,^_d) stand for the spatial 
coordinates and the microscopic velocity, respectively. The right hand side Q{f,f) is the 



collision term describing the interaction between particles. In this paper, we are focus- 
ing on the transportation part, thus the coUisionless Boltzmann equation with vanished 
Q{f, f) is considered. 

The moment method proposed by Grad [7] approximates the distribution function by 
a finite set of moments. To achieve this, we expand / into the Hermite series as in [2]: 

/(t, ^, = E /"(*' ^)^e(t,-),« i^nf^] ' (2.2) 

where a = (ai, • • • , ao) is a D-dimensional multi-index, and the basis functions are defined 
as 

'HBA^) = X{^7^^^^He^,{^d)e^v[-^-^y z = (zi,---,z^)GM^, (2.3) 

where He^ is the fe-th degree Hermite polynomial: 

i/efc(x) = (-l)^exp(^)^exp(-^), fc G N. (2.4) 

In ()2.2p . u{t, x) = (tii(t, a;), • • • , uoit, x)) and 0(t, a;) denote the macroscopic velocity and 
temperature, respectively, and they are related to / by 

p{t,x)= I f{t,x,$)dt 

Jro 

p{t,x)u{t,x)= f U{t,x,^)dt (2.5) 

p{t,x)\u{t,x)\^ + Dp{t,x)e{t,x)= I \^\^f{t,x,i)d^, 

Jrd 

where p stands for the density of the gas. The following relations can be deduced from 
the orthogonality of Hermite polynomials: 

D 

f0 = p, fe,=0, ^/2e, =0, j = !,•••, A (2.6) 

d=l 
where ej is the Z?-dimensional multi-index with its j-th component to be the only nonzero 
one and equals to 1. 

The moment system has been deduced in [5], and here we directly present the result 
therein: 

[ dt ^^ dt ^"-'^ ^2dt^ ^"^'^'^ i 

\ d=i (1=1 / 

D 






D D 



j=l d=l ^ 

I ^ ^ de 

+ 9 E E 7^ (^/a-2e,-e, + «,/a-2e, + {^j + l)/a-2e,+e,) =0, Va G N^. 






(2.7) 



In this equation, /^ is taken as zero if any components of /3 is negative. Some special 
choices of a lead to the classic hydrodynamic equations: 



where p^+e is the pressure tensoio and qj is the heat flux. They are defined as 

Pe,+ej = te - Ui){S,j - Uj)fd$, = 5ijpe + (1 + 5ij)fe,+ep (2.9a) 

J9.0 

1 /■ ^ 

g.- = 2 / 1^ - "l'(€^- - ^^O/ d^ = 2/3e, + Y. /-.+2e<^' (2.9b) 

where 5 is Kronecker's delta symbol. We refer the readers to [5] for the detailed derivation 

of (ESD. 

Since ()2.7|) forms an infinite set of moment equations which are not suitable for practical 
use, the moment closure is in need. The simplest way is to select an integer M ^ 3 and 
force /c = if |a| > Af , and the result is the Gr ad- type system with ( ^ ) moments. 

2.2 Regularization with ID velocity space 

It is well known that the lack of global hyperbolicity is one of the major defects of Grad's 
moment equations. For the thirteen moment case, the hyperbolicity region has been 
analytically obtained in [11]. The construction of globally hyperbolic moment systems 
is very meaningful to the robustness of fluid simulation using moment approximation. 
In this direction, a general method by Levermore in [^ on the construction of symmetric 
hyperbolic moment systems is proposed. Later, Torillhon [T3] raises a clever idea to enlarge 
the hyperbolicity region of the 13-moment system by using Pearson-IV-distributions. In 
[1], we have studied the general ID moment systems and found a way to make globally 
hyperbolic regularization based on the characteristic speed correction. Here we are going 
to give a brief review on the results therein. 

When D = 1, the multi- index a becomes a natural number. Substituting (|2.8bp and 
(j2.8cp into (j2.7p . we can eliminate the time derivative of the velocity and temperature. 
Thus the M-th order Grad's moment system can be written in the form of a quasi-linear 

system 

dw dw 

^ + A(«,)-=0, (2.10) 

where A is a matrix dependent on w, and 

w = {p,u,e,h,--- Jm). (2.11) 

In [1], we have obtained the following results: 



^In some literatures, the pressure tensor is denoted as ptj, i,j = 1, ■ ■ • ,D. Here the special subscript 
Ci + Cj is used to match the form of general moments fa+e. for convenience in later use. 



1. The characteristic polynomial of A{w) is 

|AI-A| = 9 2 HcM+i 
(M + 1)! 



[[^X-uf-e)fM-l + 2{\-u)fM] 



(2.12) 



2p 
2. By adding the regularization term based on characteristic speed correction 

to the right hand side of the last equation of (|2.10p . the system is turned to be 
globally hyperbolic and the eigenvalues of the regularized moment system are 

u + C,, M+i\^, j = 1, • • • , M + 1, (2.14) 

where Cj^k is the j'-th root of the Hermite polynomial Hek{x), noticing that Hek{x), k G 
N has k different zeros, which read Ci^^, • • • , C^^k and satisfy Ci^^ < • • • < C^^k- 

The second result gives a practical implementation of a globally hyperbolic regularization. 

2.3 Reformulation of the moment system 

In order to facilitate the studying of the moment system when D >2, we rewrite (|2.7p in 
another form. Let p = p9, then p = jj '}2d=iP'i(^di ^'^'^ ^^ have 

dxj p dxj Dp ^ dxj 

By substituting (j2.8|) and (|2.15|) into (j2.7p . the following equation is obtained with some 
simplification: 

^/« I V (n ^f--^ I ,, dU . f^, ^. dfa+eA ^( 9^j)\ Op 

+ 12J2^( ^fc'-e.-e, + ("i + l)/a-e,+e, " JJlPe^+e^ ) (2.16) 

+ VV \ ( -l^z^\ ^%±fd , ^^P^ 1 , ( Ca\^ dqj _ 
f^,h\^ P ) ^^^- '^Dpdx.j^y Dpjf^^dx, ' 

where Ca and Cg ^ are defined as 

D 



C„ = X]/«-2e,, (2.17a) 

k=l 
D 

C9.I = E (^/«-2e.-e, + («, + l)fa-2e,+e,) ■ (2.17b) 



fe=l 



Then collecting (^, (fZTHI) and (l2:9a)l . we get 

'^ ^ (2.18) 

+ X;i:(25.+l)/2e.-e.+e,|^+E(2c^. + l)%^ = 0, ^ = 1,..,Z?. 

The (|2.8p together with ()2.18p and (|2.16p form a moment system with infinite number of 
equations, which is equivalent to (|2.7p . 



3 System in ID Spatial Space 

In order to derive the regularization term to achieve the hyperbolicity of the moment 
systems as in Section [Z3t we first consider in this section the special case with homogeneous 
dependence of the distribution function on spatial coordinate x except for xi direction. 
Since the velocity space is multi-dimensional, the result in this section is essential different 
from [1] . The general case in multi-dimensional spatial space is studied in the next section 
based on the results herein and the Galilean invariance of the regularization. 
In ID spatial space, the distribution function f{t,xi,^) satisfies 

^7 + 61^=0, teR+, XI GM, ^GM^. (3.1) 

ot ox I 

The moment system in Section 12.31 degenerates to a simpler form. The conservation of 
mass, momentum and energy (|2.8p turn into 



dp dp dui 

D 



The moment equations (|2.16p become 

dt dxi dxi dxi 2p ^'"dxi 

y^ I fa-eg dpe^+eg ^6, a dp2ea \ _ C^dqi_ ^^ 

^ \ p dx\ 2Dp dxi j Dpdxi 

where Ca and Cg I are defined in (|2.17p . The governing equations of p2ei (|2.18p turn into: 
for i = 1, • • • , D, 

dp2ei/2 dp2ej2 , 4, . . .5^1 



dt dxi 2 dxi 

dud 

OXi OXi 



+ J2i26a + l)/2e,-e,+e.|^ + {2Sa + 1)^^ = 0. 



6 



Analogously to the moment system in Section [2. H let fa = 0, \a\ > M for M > 3, and 
then (j3.2p and (jS.Sp . together with ()3.4p form a closed moment system corresponding to 
dSl]). 

To facilitate the reading below in studying the moment system, some notations are 
introduced as follows: 



if a = (ai, • • • ,a.„) E M", then a{i:j) = (oj, • • • ,aj), 

if A = {aij)nxn G K"""", then A(i, j:A;) = (ajj,--- ,ai,fc), 



A(i:Z, j:A;) 



/ Ojj Oij + l 



V «/,. 



a«,7+i 



a/,fc J 



(3.5a) 
(3.5b) 

(3.5c) 



P = {1,2,--- ,D}, 



a = (0, a2,--- jOd), 



a 



a! 



D 

n 



D 



I 

1-) 



a,;! 



a 



Z]«*' 



i=l 



Sdm = {a G N^ I |a| < M}, Sd^M) = {/3 G 5d,m | /3 = a}. 

We permute the elements of Sd,m by lexicographic order. Then for any a G Sd,m 

holds, where Moict) is the ordinal number of a in Sd,ai, and the cardinal number of set 
'Sd,m is N = MoiMeD) = ( d ) > which is total number of moments if a truncation with 
\a\ < M is introduced. In addition, it is clear that for each a,/3 £ So,m and a 7^ /3, 



(3.6a) 
(3.6b) 
(3.6c) 

(3.7) 



5D,M(a)P|5z),M(/3) = 0, Sd,m = U SD,M{a), 



simultaneously hold. 



(3.8) 



3.1 Structure of coefficient matrix 

Similar to the ID case, a truncation with |a| < M, M > 3 is applied. Let w G 
each i,j E D, and i 7^ j, 






and for 

(3.9a) 
(3.9b) 
(3.9c) 

(3.10) 



where Am depends on ([32]), dMj) and (j33|) . 

Clearly, all the matrix Am for any D € N"*", 3 < M € N are well-defined though quite 
complex. Here we first give some simple examples and conclude a few basic properties of 
the matrix Am- 



wi = p, 

P2e, 
WAfD(2e,) - —, 

Wj\fo{a) = fa, 3 < \a\ < M. 

Combining (|3.2p with (j3.4p and (j3.3p . we obtain 

dw ^ dw 

— - + AAfT^— =0, 
ot ox I 



Example 1. If D = 2, the ordinal number of a in Sd,m is Md{oi) 



(ai+Q2+l)(Qi+a2) 



+ 



02 + 1- The permutation of w is showed in Fig. 2(a). As the simplest case, the matrix 
A3 is 



Ul 







2p 

3e/i.i 

-¥' 

2p 



whereas pi = p2ei , fn 



P 

Ul 



2/1,1 






Ul 



Pi 

/l,l 



3/3,0 
2/2,1 

/l,2 






Ul 






Ul 





2p 

2/2,0 
P 

2p 





pS-f2,0 
P 

p 

/2,0 



Ul 

2/2,0 

p 
3/1,1 

2p 



P2 
4/3,0 
3/2,1 
2/1,2 

/o,3 

= /mei+nea ■ U ^ > 3, /or any a G 



/m 

2p 






3 



Ul 





1^2 







2 



Ul 












1 




Ul 















Ul 



(3.11) 



and 3 < |a| < M, 



Aa/(1:10, 1:10) = A3, 
AA/(7Vz)(a),7VD(a)) = ui, 
AM{J\fD{a),J\fD{a-ei)) = e, ifai>0, 
AA/(AAD(a),A/'D(a + ei)) = ai + 1, if \a\ < M, 

AMiMoia), 1:9) = (-^<^S, ^/a-2e, + (ai + l)fa - ^P2e„ 



0fa 



+ (ai + l)/< 



a+ei— 62 



2l2z^ + 



a 



(1) 



/, 



a— 62 



a 



2^^61+62, 

(1) 



(3.12a) 
(3.12b) 
(3.12c) 
(3.12d) 



(3.12e) 



3Cq, 



p 2p p 2p 2/9 

where Ca and Cg ^ are defined in equation (|2.17p . We remark that 
• an entry Aj\,/(z, j), if not defined above, is taken as zero; 



, 0, 






), 



for \a\ = 4, some Ai[j{i,j) may be double defined in (j3.12cp and (j3.12ep . the value 
of which is the sum of the both expressions. 



Fig. [T] gives the sparse matrix pattern of Am with M = 8 and D = 2. It is observed 
that there is no more than one nonzero component of AA/(i,i + 1 : N), for each i = 1, 
• • • , N. Precisely, there is a unique nonzero component as 1 < i < MdUM — l)eci) 
and AM{i,i + 1 : N) = as i > MdUM — 1)6/)). Noticing that the column index of 
the nonzero entries in the upper triangular part of Aj\,/ on different rows are different 
from each other, this makes one recall the form of lower Hessenberg matrix, of which the 
only nonzero entries in the upper triangular part on the i-th. row is located at position 
{i,i + l). This property of Hessenberg matrix makes it very convenient for one to calculate 







10 



column 
20 30 



40 



-AAf(l:10, 1:10) 



AM{J^D{ct),AfD{a + ei)) 

Am{Nd{cc),Md{c,)) 

AM(A/'D(a),A/'D(a-ei)) 




Figure 1: The sparse matrix pattern of Am with M = 8,-D = 2. Its nonzero entries are 
defined as in ()3.12p . 



its eigenvectors once eigenvalues are given using a row by row sequential procedure. Here 
Aj\// is essentially the same as lower Hessenberg matrix on this point, and we notice that 
its lower triangular part is sparse, hence we are provided the approach to calculate the 
eigenvalues Ajv/ together with the corresponding eigenvectors using the same technique. 

Furthermore, (|3.12|) shows the diagonal entries of the matrix Ajvf are all ui, and the 
entries of the matrix A^ — "Uil are independent of Uj, i £ D, where I is the N x N identity 
matrix. In fact, (j3.10p can be written as 



— + (Am - ^il) ^ 



0, 



(3.13) 



D 



where — is material derivative defined as 
Dt 



D 
Dt 



d_ 
dt 



ui 



d 

dxi 



Hence, that Am — uil are independent of Uj,i S P indicates the moment system is 
translation invariant. On the other hand, eigenvalues of Am can be written in the form 
ui + a, where a is indeterminate and independent of Ui, i & V, and eigenvectors of it is 
independent of Uj, i £ V, too. 
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(a) The permutation of ly 



(b) A permutation of lu' defined in example [3] 



Figure 2: The permutation of moments while L) = 2, M = 8. Each node stands for one 
moment. The marks in the lower right of the node shows the expression of the moment, 
while the number in the upper left represents the ordinal number in w or w' . The dashed 
arrows depict the path of the corresponding permutation. The left one is the permutation 
of to, and the right one is a permutation of w' defined in example [3j 

Example 2. Considering the case D = 2, we can write out the matrix Am according to 
example\^ for any 3 < M G N. If we let /« = for all a G N^, and \a\ < M except /o 
and fMen direct calculation gives the characteristic polynomial of Am as 



lAI- A 



M\ 



X HCM 



X — ui 



X — ui 



ri/2 



Ve 



e^l^ + {-lf"-'M\fMe. 



(3.14) 



X HCM+I 



ihr) ^^''^'^" + (-1)''"'(^ + mMeA^ - "l) 



The matrix Am, obviously, has complex eigenvalues, for some fMei- 

Analogously, involved calculations with the help of computer algebraic system show 
that the matrix Am has complex eigenvalues for some admissible w, for any D >3. This 
reveals that Am is not diagonalisable with real eigenvalues for some w as D = 1 in [1]. 

Remark 1. Since moments fa are related to f{t,x,$,) by ()2.2p . the moments fa can not 
be arbitrary number if the distribution function is to be kept positive. Particularly, p and 
9 given by ()2.5p clearly satisfy 

p>0, 9>0. (3.15) 

Though (|3.15|) is not enough to provide us a positive /(t, x, ^), the discussion in this paper 
requires no further constraints on the other moments. Hence, in this paper the admissible 
w stands for w which satisfies ()3.15p . 

Example 3. Actually, it can be observed that the matrix Am is reducible if we rearrange 
w as w' using another permutation rule. In case of D = 2, the rule reads: 



10 



1. The moments with 02 < 2 are arranged at first using the lexicographic order; 



2. The rest moments are arranged then using the lexicographic order based on index 
transformed as (a2,ai). 

Clearly, w and w' are related by a permutation matrix P that w' = Pw. The Fig. 
\2(b)\ gives a schematic diagram of the permutation rule for w' with M = 8. Let A^^ = 
PAAfP-\ then 

^^' " ^^'-^ (3.16) 



+ A 



0, 



dt ' "^dxi 

holds. Fig. \3.1\ gives the sparse matrix pattern 0/ A^j. It is obvious that A^j is reducible 
(see e.g. fEl for the definition), and can be reduced into M—1 blocks. 5D,A/(ei) U SD,M{e2) 
U •Sd,m{'2(^2) is one of the blocks, and S^j{a), for each a G N and a / ei, 62, 2e2 is another 
block. 
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Figure 3: The sparse matrix pattern of A^j with A-I = 8,D = 2. A^j is reducible. 

These examples show some very useful properties of the matrix Am as follows: 

Property 1. Matrix Am satisfies the following properties: 

U](l). For each a G N^ , \a\ < M , let i = Moioi), then there are no more than one entries 
ofAM{i, i + l:N) nonzero. In fact, there is a unique nonzero component as \a\ < M 
and Ajv/(^, i + l:iV) = as \a\ = M . 

\^2). The diagonal entries of the matrix Am are allui, and entries of the matrix Am — uil 
are independent of Ui,i G T>. 

n](3)- AA/(if ) may be not diagonalisable with real eigenvalues for some admissible w. 
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U}(4)- ^M is reducible, and can be reduced into ( £)_i ) — 2(D — 1) blocks. 5D,A/(ei) IJ 

Q SD,M{ek) I U ( U SD,M{^ek) ) is one of the blocks, and Sd,m{oi.), for each 
\k=2 / \k=2 / 

a G N^""*^ and a ^ ei, e^, 2efc, A; = 2, • • • ,D is one of the blocks. 

3.2 Globally hyperbolic regularization 

For ID case, the regularization as ()2.13p was proposed such that the moment system turns 
out to be globally hyperbolic (see [T] for details). Actually, the regularization therein can 
be extend to multiple dimensional systems. For D G N'^, let us start from the definition 
as below: 

Definition 1. For any \a\ = M , let 

D 

where 

d=l ^ \d=l / ^ 

As in ID case f]}/, 7^a-/,d(ck) is the regularization terms based on the characteristic speed 
correction. 

With (J2J5|) . we have that 

For the case that the dependence of / on a; is only on xi, we have that 

^m,d(«)=0' for j = 2,--- ,D. 

This leads to TZm,d{oi) = TZ\,jjj{a). The regularized system is obtained by subtracting 
'T^m,d{oi) from the governing equation of fa in ()3.10p . for \a\ = M. 

Definition 2. Aj\/ is called the regularized matrix of the matrix Am, if it satisfies that 
for any admissible w, 

J dw dw sr-^ 1 / \ 

d^ ^ dxT ~ ^ i^MMoinMoia)^ (3-19) 

-*- -*- \a\=M 

where I^ is the k-th column of the N x N identity matrix. 

The regularization terms only change a few entries of the lower triangular part of Aj\/ , 
with the order of the corresponding moments equals to M, so that the Properties \T[1)\ 
[112)1 WlA)\ of Am are also valid for Am, while the Property \T\3)\ is changed to be the 
diagonalisability of Am over the real field M. Actually, we have the following theorem: 
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Theorem 1. The regularized moment system. 

dw 2 dw 

is hyperbolic for any admissible w. 

The definition of the hyperbolicity shows that this theorem is equivalent to the diago- 
nalisabihty of P^m with real eigenvalues for any admissible w. Before proving this result, 
we first make some simplifications and give several useful lemmas. 

Let us denote that 

d = {dj)Nxi, di = p~\ (3.21a) 

dj+i=d-^/\ j = V, ^AT^M = /5"'^"'"'/', 2<\a\<M, (3.21b) 

A = diag{di,d2,- ■■ ,dN} , (3.21c) 

p,,+,, = ^, i,k = V, 9a = ^\. (3.21d) 

By the virtue of Property ^2}] we let 

Am = ^il + VoA-^AmA. (3.22) 

Then, we can obtain properties of Am by studying the matrix Am- Since Am is related 
to Am — uil by a similarity transformation by a diagonal matrix, Property flTl)! holds 
for Am- Hence, it is convenient to calculate the eigenvectors of matrix Am- Firstly, we 
denote by some symbols as 

if ai = 0, (3.23) 

rp,^,/2 = yr„ (3.24) 

if at least one aj < 0,j £ V (3.25) 

k £ V\{1}, (3.26) 



^'"jVo(c-) 


-mD-i{&)^ 


Tm 


= >^rp, 


"■fn 


= ru = 0' 


^P-i+-k 


= ^rui^ , 


rh.. 


d=l 


'^^AfDic) 


-'^{ru^Oi. 



(3.27) 

(r/^ + G{a)) - G{a), if ai / 0, |a| > 3. (3.28) 

where fA/'o_i(a) ^^^ A are indeterminate parameters, and 

G(a) = J2 5-e,r„, + f E ^ - ^) E 9a-2e, - (3-29) 

d=l \d=l / fc=l 

For better readability, here we adopt the notations [e.g. Wk, fa, f2ei] as the subscript 
of r. This does not mean the subscripts are taken as the value of them, but only taken 
literally as the notations themselves. We collect these r^^j,, with Wk as a component of w 
and k = 1, - - - , A^, to produce a vector r G M as 

where N is total number of moments. Here it is clear that r is prescribed once f_/v^_j(^) 
and A are all given. With particular setup of these parameters, A and r is turned out to 
be a pair of eigenvalue and eigenvector of Am- Precisely, we have the following lemma: 
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Lemma 1. r j^ is the right eigenvector of the matrix Am for the eigenvalue A if 

%l'tiy i'-^+G{a))=0 (3.31) 

holds, for all \a\ = M . 

Proof. Let i = A//) (a), with \a\ < M, then we need only to verify 

AMii,l-N)-r = \r^^ (3.32) 

always valid. Since Am is determined by ()3.2p . (|3.3p and (|3.4p . and Am and Aa/ are 
defined as (|3.19p and (|3.22p . respectively, we can write any entries of Am- Now let us 
verify the equation (|3.32p case by case: 

• For a = 0, 

AM{i,l-N)-r = l-ru,=Xrp. (3.33) 

• For a = ei, 

AM{i, l:N)-r = 2- r^^^j^ = AV^ = Ar„,. (3.34) 

• For a = Ck, k = 2, ■ ■ ■ , D, 

AM{i, l:N)-r = l- r^^^^^j^ = Ar„,. (3.35) 

• For a = 2ei, 

3 ^ 

AM(i, l:N)-r= -r„^ + ^ 'igo.+e^.e^ru^ + 2,rf,^^^ 

d=l 



^'^wi +%2ei?'«i +3 f ^ rp-g2e^r, 1 *- • '' 

A^ 



• For a = 2ek, k = 2,- ■ ■ , D 
AM{i,l-N)-r 






- ^^P2.,/2- 

For a = ei + efc, A; = 2, • • • ,D 
AM{i,l-N)-r 



1 ^ 



d=l 



1 Ar _1 ^ _ _ 
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(3.37) 



(3.38) 



For 3 < |a| < M, ai > 



1 -(1), 



AM(i, 1 : iV) • r = 1 • r/_^^ + {a, + l)r^„^^^ - -C^;>, 






1 ■ 



= ^1+^2, 



whereas 

D 



ail 



rj^ + G{a)) - G{a - ei) - (ai + l)G(a + ei). 



Ja~2ek, 



(3.39) 



k=l 
D 

^e,a = 2^ {9a-2ek-ei + ("l + l)9Q-2efe+ei) , 
fc=l 

D D 

fe=2 fc=l 

Xi = r/^_^^ + (ai + l)r/^^^^ , X2 is the rest terms. 
Substituting ()3.29p into Xi yields 

+ (ai + 1) (^^^ (^/a + Gia)) - G{a + ei)) (3.40) 



For X2, using (|3.'26p . we get 

D 

^2 = ^ 2^ (5a-ei-2efc + ("l + l)ga+ei-2ek " Aga_2eJ 

fc=l 
D 

+ X] '^"rf i9a-ei-ea + ("1 + 1)50+61 -e^ " Afifa-eJ (3.41) 

(D ^ \ D 

Yl '^Jj~ Yl (5'"-ei-2efc + («1 + l)5'a+ei-2efc " >^9o 
d=l J k=l 
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Now we calculate G{a — ei) + (ai + l)G{a + ei) — XG{a). Some simplification gives 
G{a - ei) + (ai + l)G{a + ei) - XG{a) 

D 



d=l 

/ ^ r \ ^ 

+ Yl '^ ^ X (5a-ei-2efe + ("l + 1)5^+61-26^ - ^9c,-2ey,] 

\d=l J k=l 

(I330|), (fOT]) and ([3:12]) show 

Xi + X2 = A:^^^(r;, + G(d)) - \G{a) = Xr^^. 



(3.42) 



(3.43) 



For 3 < |a| < Ad, ai = or a = e^ + Cj, j > k > 1, it is the case that to let 
''/a-e — ^ ^^^ (|3.39|) . which is actually part of p.25p . Hence, (j3.32p is valid in this 
case. 

For \a\ = M, if ai > 0, then (j3.19p and (|3.22p show this case equals to let 

^/.+ei + Yl 9a-ea+e^ru^ + Yl 9a-2e^+e^ Yl jfv^nl'^ ~ t" A " ^ ^^'^^^ 

d=\ \d=\ J \%=\ ) 

in (fOQD . Since 

D / D \ / D 



P ' 



G(a + ei) = Yd'^-^d+^i^^d + X]^"-26,+6i X] ;d^P2^'/^ ~ 2 

(i=l \d=l / \J=1 / 

we need only to prove 

Actually, it is what (|3.3ip tells. 

If «! = 0, (|3.19p and (|3.22p show this case equals to let rf^_^ = and (j3.44p valid in 
()3.39p . The former is part of ()3.25p . while the latter is proved above. Hence, ()3.32p 
is valid in this case. 

Collecting all the cases above, we conclude (j3.32p is valid for arbitrary a. The lemma is 
proved. D 

For any a, let /3 = a + kei,k E N, then a = /3 holds. Therefore, r^-, \a\ < M is 
equivalent to r^^, \a\ = M. Hence, in the lemmalU parameters r^, |a| = M and A, are all 
indeterminate. Let A^^ = MD-iiMeo), and 



V = {vi,V2,--- ,VnJ G 



t>N^ 



Since r is determined by v and A, by studying the space of the parameters v and A, we 
can fully clarify the structure of the eigenvectors of Am- We have the following lemma 
that 
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Lemma 2. Am has N linearly independent eigenvectors. 
Proof. For |a| = M, ()3.3ip can be written as 

• if d = 0, then A/'d_i(q;) = 1 and 

viHeM+i{y) = 0, 

• if a = e^, k £ P\{1}, then A/z)_i(a) = k and 

VkHeuW = 0, 

• if a = 2ek, k £ V\{1}, then 



(3.45) 



(3.46) 



D 



0= {^^n-ii2e,)+Yl 



""-^^-^-^UeM-i(A) 



D 



VAfD-i{2&k) +Y1 



D 



VAfo^ri2ea) A^ n\ rr /xN 



(3.47) 



\ d=2 

• if a = ek + ei, k j^ I, and fe, / E P\{1}, then 

^Aro-i{Q)^eA/_i(A) = 0, 

• otherwise (3 < \a\ < M), 



(3.48) 



D 



D 



D 



VMd 



-l(") + Y^ da-e^Vd + ^ 5a-2e, X] 



WD-i{2ed) A^ ■^l n W (W n 



d=2 



\d=2 



D 2D 



Let 



2A = (^eM+i(A), ••• , i?efc(A), ••• , /^efc(A), •••,Fei(A)), 



(3.49) 



(3.50) 



(''m+Y-\t') entries 



#{a I |a| = M, ai = A: — 1} 



where the A/'D-i(a)-th component of z;^ is ii^ea^+i(A), and the cardinal number of set 

'D-l + M -k^ 
M + l-k 

Equations ([335]), (laliHjl . (13:17)1 . (IHIiHjl and (I3:i9]l can be collected as 

ZAoBti = 0, (3.51) 

where c = a o b, stands for Cj = Ojftj, i = 1, • • • , n, and B is a {N^ + 1) x (N^ + 1) real 
matrix. Precisely, the formation of B is as 



B 



10 

B21 B22 

B31 B32 I 

t t t 



D row 

D(D~l)/2 row 
JV„-D(Z)+l)/2 row 



(3.52) 



D col £i§^ col 7V„-£i§+ii col 
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where I is identity matrix, whose dimension is context depended. The first D rows of B 
are arising from (|3.45p and (j3.46p . the following D[D — l)/2 rows are arising from (j3.47p 
and ()3.48p . and the rest Ny — D{D + l)/2 rows are arising from ()3.49p . 

The properties of B are here further clarified. We denote the entry of B located at 
(i, j) position as hij, whereas i,j =0, • • • , N^. Noticing that entries of B21 are from ()3.47p 
and ()3.48p . we have 

A^ 1 



0, 



else. 



And (|3.47p and ()3.48p make that B22 is a summation as 



B 



22 



i + j^a 



(3.54) 



where fi G 



D(D-1)„D(D-1) 



is a matrix with its entries ujij as 



UJ. 



V 



1, if i = A/'D-i(2efc) - D, j = No-ii'^h) - D for some k, I G P\{1}, 
0, else. 



Since there are at most D — 1 nonzero entries in each row of $7, it is clear that B22 is 
strictly diagonally dominant thus nonsingular, and then B is nonsingular. With the entry 
value of fi as given above, one can check fi^ = {D — l)d holds. Hence we can get the 
inverse of B22, which reads: 



b: 



-1 



1 



-Cl. 



(3.55) 



22 2D-1' 

Meanwhile, since B is a nonsingular block lower triangular matrix, we can get its inverse 



as 



B 







B22 B21 



B 



22 





I 



D row 

D(D-l)/2 row 
N^-D(D+l)/2 row 



(3.56) 



Let 



t t t 

D col 2i^ col N^^^i^ col 
B = diag{I, B22, I}, 



(3.57) 
be the diagonal blocks of B. The inverse of B is B~^ = diag{I, B22 , I}, and we have 

BB~' = -B^i I . (3.58) 



(3.59) 



I 








B21 


I 





* 


* 


I 



Since B is nonsingular, for an arbitrary j € {!,••• , Ny}, we let 



^j' 



where Ij is the j'-th column of the N^ x N^ identity matrix. Actually, v'^^' is the j'-th 
column of BB~^. Notice in (j3.50p for any a that \a\ = M, A/£)_i(a)-th component of z\ 
is i7eQ,j-|_i(A). For the a satisfying 



a 



M and j = Md~-i{ci), k = ai + 1, 
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we choose A such that 

HekiX) = 0. 

Then we have that the j-th component of zx vanishes 

zxj = Hek{X) = 0. 

This makes ()3.5ip vahd that 

BB-ii;(^) = plj, ZAj =0, j = 1, • • • , iV.. (3.60) 

Since r is depended only on v and A, we denote r^, j to be the vector prescribed by 
the given v^^', j = A/z)-i(a), and A = Ci^k, when k = ai + 1, for arbitrary |a| = M, i = I, 
■ ■ ■ , k. It is clear that Ci^k and Tq^j are a pair of eigenvalue and eigenvector of Aj\/ that 

The eigenvectors of Aj\,/ can be divided into a cluster of classes, each of which is: for 
arbitrary |q;| = M, 

{r^^i I i = I,... ,k, k = ai + l}. 

This fact essentially stems from the reducibility of the matrix A^- 
Notice that 

1. The components of v are a subset of r's components, linearly independent v^^'^s 
determine linearly independent r's; 

2. Eigenvectors belongs to different eigenvalues are orthogonal and the k zeros of Her- 
mite polynomial Hek{\) are different. 

We have that r^ j, i = 1, • • • , /c, when \a.\ = M and k = ai + 1 are linearly independent 
and the matrix Am has 

v^^ fD-l + M - k\ fM + D\ , , 

linearly independent eigenvectors. 

On the other hand, for arbitrary ra,i, there exists a unique /3 satisfying /5 = (i— l)ei +d, 
hence, there is a one-one mapping between r^ j and a with \a.\ < M. So we can also get 
Aj\,/ has N linearly independent eigenvectors. This completes the proof. D 

With the help of Lemma [21 it is not difficult to get the following result: 

Lemma 3. Let 

Vi,m = Hem.+i{X), meN, (3.62) 

m 

VD,m = n ^D~i,k, K £> G N+. (3.63) 

fc=0 

'Pd,m is the characteristic polynomial of Am- 
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Proof. In case of D = 1, the result has been proved in [T]. Here we give the proof for 
D > 2. In the proof of lemma [2l (|3.60p shows that the characteristic polynomial of Am is 

A/+1 M+1 

Vd,m= n ^efe(A)i i^^+i-^ ) = II HekiX)^ 0-2 ). (3.64) 

fc=i fc=i 

Now we need only to prove that Vd,m satisfies ()3.63p . Here we use induction argument 
on D. As D = 2, ()3.64p can be written as 

M+l M 

V2M = n H^kW = n A,m. (3.65) 

fc=l m=0 

We assume that (|3.63p holds for D — 1,D > 3. With the induction hypothesis, we have 

(D-\-l + m-k\ \ 
U D-1-2 I 



M M /m + l 

n ^D-i,^ = n n ^^^w 

m=0 m=0 \fc=l 

M+l 

= n Hek{\) 



M+l 

Y^M jD-2+m-k\ 



(3.66) 



' D-l + M-k\ 
D-2 ) 



M+l 
= n Hek{\) 
k=l 

= 'Pd,m- 
This completes the proof. D 

With the relation of Am and Am (|3.22p . we have the following theorem. 
Theorem 2. Let 

V,,,n = Hem+1 (^^) ^^"''"'^^'' '^ e N, (3.67) 

m 

rD,m = n T'D-i,k, l<Den+. (3.68) 

fc=0 

'Pd,m is the characteristic polynomial of Am- And Am has N linearly independent eigen- 
vectors, which read 

Ta^i = A~ ra,i, for eigenvalue Xi^k = ui + Cj^feVP (3.69) 

for all \a\ = M , i = 1,- ■ ■ ,k, whereas k = ai + 1. 

Proof. Since 



^M = Ui 



I + VoA-^AmA, (3.70) 

and A is nonsingular, so any r^^j G M^ is the eigenvector of Am for the eigenvalue Cj^^, 
then r^^i = A~ r^^i is the eigenvector of Am for the eigenvalue ui + Cj^fcV^. Using 
Lemma [1] and discussion in Lemma El we obtain p.69p . Lemma [2] shows Am has N 
linearly independent eigenvectors, so A^j also has N linearly independent eigenvectors 
and ()3.69p gives a set of basis. 
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Lemma [2] and ()3.69p show that the characteristic polynomial of Aj\/ is 

M+l . . , ^ X (D-l + M-k\ 

Similar as that in the proof of Lemma [3l Vd,m is thus the characteristic polynomial of 
Am- □ 

Theorem [T] is now straightforward: 

Proof of Theorem [IJ With Theorem [2l we declare that Am is diagonalisable with real 
eigenvalues directly, that is, the moment system ()3.20p is hyperbolic. D 

4 System in Multi-dimensional Spatial Space 

As the main result of this paper, here we give the general hyperbolic moment system 
containing all moments with orders lower than M. Without the assumption that the 
dependence of / on X2, • • • , xd is homogeneous, according to the discussions in Section 
12.31 Grad's moment system can be written in the following form: 

dw v^ , ^ , , dw ^ , ^ ^ , 



where w remains the same definition as the one-dimensional case ()3.9p . and Mj, j = 
1, • • • ,D are square matrices depending on w. Comparing with (j3.10p . one immediately 
has Ml = Aj\/. Similar as Definition [21 we give the following definition: 

Definition 3. For j = 1, • • • ,D, Mj is called the regularized matrix of the matrix Mj, if 
it satisfies that for any admissible w, 

^fi/7 /'tin 

^^d^ = ^'d^~ ^ ^m,d(«Ka^z,H' (4-2) 

^ ^ \a\=M 

where Ik is the k-th column of the N x N identity matrix. 

Now the multi-dimensional regularized moment equations can be written as 

ljf + E'^.Wa^ = o. (4.3) 

Recalling that 

D 

7^M,D(a) = J^7ei,B(a), (4.4) 

i=i 

one finds that the multi-dimensional regularized moment system is obtained by subtract- 
ing TZM,D{ci) from ()4.2p for all \a.\ = M. Applying such an operation on ()2.7p . we can 
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reformulate the regularized moment system as 

[ dt ^^ dt ^"-^^ + 2 at ^ •^""'^'^ i 

V (1=1 d=l / 

D 



+ Zl Zl 7W^ {^f»-^d-<^j + ^i/a-ed + (1 - '^|a|,A/)(ai + l)/a-ed+e,) 



1 ^ ^ 5^ 
+ 2 2]] Z] 5^^ (^/a-2erf-e, + %/a-2erf + (1 - '5|a|,A/)(aj + l)/a-2erf+e,) = 0, 

|a| ^ M. 

(4.5) 

Actually, ()4.5p is away from (j4.3p only by a linear transformation due to ()2.8p to eliminate 
the time derivatives of u^ and 0. Precisely speaking, there exists an invertible matrix 
T^{w) depending on w such that (|4.5p is identical to the following system: 

T(«^)^ + EtMM,(^)|^ = 0, (4.6) 

If we let all partial derivatives with respect to Xj with j > 1 to be zero, ()4.5p reduces to 
the one-dimensional hyperbolic moment system (|3.2Up in Section O Comparison of (|4.3p 
and ()3.20p clearly shows that Mi = Km- 

The following theorem declares the hyperbolicit}o of the multi-dimensional regularized 
moment system ()4.3p : 



Theorem 3. The regularized moment system ()4.3p is hyperbolic for any admissible w. 
Precisely, for a given unit vector n = (ni,--- ,n£)), there exists a constant matrix R 
partially depending on n that 

D 

^njMj(w) = R"^Aa/(Rw)R, (4.7) 

and this matrix is diagonalizable with eigenvalues as 

u-n + Cn,mV0, l^n^m^M-Fl. (4.8) 

Actually, this theorem gives the rotation invariance of the regularized moment system 
and its globally hyperbolicity Since the translation invariance of the system is apparent, 
it is concluded that the regularized system is Galilean invariant. Precisely, if another co- 
ordinates (xi, . . . ,xd) are chosen and the vector n is along the xi-axis, then the rotated 
moment system is equivalent to the original one. This result is easy to understand: on one 
hand, Grad's moment system is rotationally invariant, since the full M-degree polynomi- 
als are used in the truncated Hermite expansion; on the other hand, our regularization 



^For multi-dimensional quasi-linear systems, we refer the readers to [8] for the definition of hyperbolicity. 
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is symmetric in every direction, which can be considered as "isotropic" in some sense. 
However, a rigorous proof of this theorem is rather tedious. 

In the hterature, two types of indices have been used in the moment methods. In 
Grad's paper [T], indices such as 

^=(^i,---,^„)GP'" (4.9) 

is used to denote the m-th order moments, while in [2], the symbols 

a = (ai,--- ,az)) gN^ (4.10) 

is used as the subscripts of |a|-th order moments. The former is convenient for mathe- 
matical proofs, while the latter is easier to use in the numerical implementation, since for 
()4.10p . the map from the index set to the moment set is a bijection, while this is not true 
for (|4.9p . If (j4.1Up and (|4.9p represent the same moment, then one has 

a = ei9i H ^e^m' m = \a\. (4-11) 

Below, both types of indices are needed in the proof of rotation invariance, and we will 
always use the variant forms of Greek letters such as •& and ip to denote the Grad-type 
indices, and normal Greek letters such as a and (3 will be used to denote indices like ()4.10p . 
The Greek letter "sigma" denotes the conversion between them. Supposing (|4.1ip holds, 
we write 

a = a{^), T? = ?(a). (4.12) 

That is, the normal form of sigma (t(-) converts indices like ()4.9p to indices like ()4.10p . 
and the variant form of sigma ?(•) does the inverse conversion. Note that for a given a, 
the Grad-type index •& satisfying ()4.1ip is not uniquely determined. Define 

1(a) = {^ G pl"l I a{'d) = a}, (4.13) 

and then in most cases, l^(a) has more than one element. For example, if D = 2 and 
a = (2,2), then 

1(a) = {(1,1,2,2), (1,2, 1,2), (1,2,2,1), (2,2, 1,1), (2, 1,2,1), (2, 1,1,2)}. (4.14) 

Thus (,{ol) has multiple values. However, there is always one special element {} E ^(a) 
satisfying 

^i^---^%|, (4.15) 

and we use this element as the value of ?(a). It is easy to find 

o-(?(a)) = a. (4.16) 

Additionally, we use o"j(t?) to denote the i-th component of cr(t?). 
Based on these symbols, we have the following lemma: 

Lemma 4. Suppose a G N and F{-) is a function on V^. If F satisfies that F{ip) is 
zero when 0"j((^) < Oj for some i £ D, then the following equality holds: 

E ^(^) = E E p(^)- (4.17) 

|/3|='fn— |a| 
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Proof. It is obvious that 

1= \J I(/3 + a)cP"^, (4.18) 

|/3|=m-|a| 

and there are no duplicate elements in the union since T{/3 + a) n T{f3 + a) = if /? 7^ /3. 
Thus it only remains to prove that ip £ I ii 

ipeV"^, and ai{if) ^ ai, Vi G P. (4.19) 

This is true since (p G T{/3 + a) for /3 = cr{(p) — a. D 

As a special case of Lemma HI we set a = and have 

E ^(^) = E E ^('^)- (4.20) 

|/3|=m 

Here -F(-) is an arbitrary function on D"^. 

Some more symbols are introduced as follows. All ?Ti-permutations of the set {1, • • • , n} 
form the following set: 

A^ = {w = (tZ7i, • • • ,-a7m) G {1, • • • ,n}™' | rui 7^ Wj if i / j}, \/m,n G N, n^ m, 

(4.21) 
which contains n!/r7i! elements. Thus when we want to construct a short vector using the 
components of a long vector, we will use the following notation: 

^^ = (^^i,--- ,^^™)e^'", V^GP", tuG^::^. (4.22) 

The remaining part is denoted as 'Q\&vj- For example, \i ■& = (1,3,2,3,1,2,1) and w = 
(5,2,4), then 

^^ = (1,3,3), i9\^^ = (l,2,2,l). (4.23) 

Below, G = {gij)DxD stands for the rotation matrix, and we suppose G is orthogonal 
and its the determinant is 1. Define 



n 



ngi^,v) = llg^,^,, v^,(^GP^ (4.24) 



i=l 



and then we have the following lemma: 

Lemma 5. For a given matrix G and multi-indices a, /3 ^ N , the following equality 
holds for arbitrary "d G "DI^I+I^I.- 

^ ^"^(^'^^ = ^ ^ n,(^^,,(/3)) Y. n,(^\^.,c^). (4.25) 

Proof. We first consider the case |/3| = 1. Suppose /3 = e^, and then ()4.25p becomes 
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For 93 G E(a + e^), one has cr{ip)\ = (a^ + l)a\. Thus ()4.26p is equivalent to 

H+i 
(a, + l) Yl n,(z?,^)= J^g^,, ^ n3(^\^„(^). (4.27) 

ipG (d+Ed) i=l <pe (a) 

For an arbitrary tp G T{a + Sd), if v^j = d, then ng('!9,c/5) = g^.dllg{'&\'&i,<f\<fi), and 
C/9\c^j S E(a). Since there are [ud + 1) choices of i such that ipi = d, the product Hgii}, ^p) 
appears (a^ + 1) times in the right hand side of (|4.27p . This proves (j4.26p . 

Suppose the lemma holds for |/3| = m — 1, and we are going to prove the case \/3\ = m. 
In order to use the technique of induction, we choose d € {1, • • • ,D} such that (3(i > 0, 
and let j3' = (3 — ed- Thus |/3'| = m — 1. Applying (|4.26p . one has 

./JS (a+/3) ^ ^ (pe {Q+/3'+ed) ^ ^ 

= "^(^y " Y. S^^d Yl ng(??\^i,(/p) (428) 

Defining i?* = ^\'&i, and using the inductive assumption, one obtains 

E ^n,(^,v.) = ^ E <7... E n,(C,^(/3')) E n,(^AC,¥') 

(4.29) 

It is evident that the right hand sides of (|4.25p and (|4.29p are the same. Thus the lemma 
is proved. D 

Now let us start the rotation. We first define 

D 

Xi = Y^gijXj, i = l,--- ,D, (4.30) 

i=i 

and denote by p, ii, 9 the density, macroscopic velocity and temperature in the new 
coordinates x = (xi, • • • , xd). If we define ^ = G^, then the orthogonality of G shows 

P= I f{Odl= I /(^)d^ = />, 

pu= f ifiO d| = / G^/(^) d^ = pGu, (4.31) 

and it follows immediately that 

D 

9 = 9, Ui = Y,9ijUj^ i = l,---,D. (4.32) 

i=i 

25 



Now we consider the general moments /q, in the coordinates x. Define z = {^ — u)/\/9 
and z = {^ — u)/\'0. Then z = Gz. The orthogonahty of Hermite polynomials gives 

(27r)^^l°l+^ 



fa 



(4.35) 



/ f{u + y/0z)ne,a{z)exjp( ^ ) dz, 

/a = ^^n / /(^ + ^^)'HeA^) exp -^ dz, 

Oil Jrd V 2 y 

From the definition of Hermite polynomials ()2.4p . it is easy to find that ()2.3p can be 
rewritten as 

?^,,„(z) = (-l)H(2vr)-Tr^|_exp(^-^j . (4.34) 

Applying the chain rule of differentiation, we obtain 

?^,,„(z) = (-l)H(2vr)"f0-^^ ^ n,(,(a),c^)^^exp(-^ 
= Yl ng(?(a),(/p)?^e,^(^)(z). 
Collecting (|T33]1 and (|05]) . one has 

-^^^"' (4.36) 

= E 3^n,(,(a),^)^(,). 

As in (|3.9|) . all the rotated moments can also be collected into a vector denoted as w. The 
equations ()4.36p and (|4.32p directly give the following result: 

Lemma 6. Based on the expressions of the rotated moments (|4.36p and (|4.32p . the fol- 
lowing equalities hold for arbitrary a G N .' 

X^ ~ dfa V^ V^ (^ivV-jr t f ■, \ ^f'^if) U'i^\ 






D „ r D 



Y^^' + ^)^— = Y L ^^n,(,(a),v.)-(Mv^) + i)^^^, 

d=l d=l tpeT>\"\ 



(4.38) 



where craif) is the d-th component of a{ip). 
Proof. Using (|4.36p and ()4.32p directly, we get 

E-t - E E ^n,M„)„) . E««,%f , (,3.) 

d=l d=l ip£j}\°'\ j=l 



Equation ()4.30p shows that 

E^^.7^- (4-40) 



d ^ d 



•^^i ^ ^^rf 
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Thus (|07)1 is the direct result of ()i:H9]) and dOO]) . 
The proof of ()4.38p is also straightforward: 



dfa+ea _V^/ , IN Y^ 0-(v?)! x^/a(¥:.) 






= Y.Y.9dj E a! "^('^(")'^) g^, (4-41) 

d=l j=l ipeV\"\ 

= E E ^n,M„).^).(.,M + i)^%^. 

This equality is identical to (j4.38p . D 

Using Lemma [5j it is not difficult to prove the following lemma: 
Lemma 7. The following equalities hold for arbitrary a G N^ : 

Ef^/"-- = E E ^n, (.(„), ^)^/.,„_.„ (4.42) 

d=l d=l(pGl)l"l 

Ei/"-- = E E ^n,(.(„),^)f /„„_,., (4.43) 

d=l d=l ip^V\"\ 

E%f -E E ^n,(.,„),.)%b^. (4.44) 

EEfF/"— '.=EE E ^n,{.(a),^)|^/„„_,_.„ (4.45) 

j:j:§j"--".-', = EE E ^n,(,(„).^)|^/„„_,.,_.,. (4.46) 

i=l d=l -^ j=l d=l <^gX)|a| -^ 

Proof. Recalling that /^ = if /3 has a negative component and using Lemma HJ we have 
the following equality: 

= E E E ^n.w°>.-)%^/^' 

d=l /3gN° ve (/3+ed) 
|/3| = |a|-l 



Let -& = <;{a) and use (|4.26|) . and we have 



^ = E E ^Es^^^ E n,W.,v.)^/,. (4.48) 

d=l i3(zfqD i=i ^g (^) 

|/3| = |q|-1 
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Now we employ ()4.20p to join two of the summation symbols in the above equation: 

^=E^i;^-^'^ E ^n,w.,^)/.(,). (4.49) 

d=l i=l (pgDl^l-i 

Using ()4.36p and ()4.32p again, we get 

D \a\ „ \a\ ^ rs~ Dr.- 

Thus ()4.42p is proved. The equation ()4.43p can be proved in a similar way. Setting 
'd = <i{a) and using Lemma HI Lemma O and (j4.36p . we obtain 

d.=l (^gX)|o! 

-E E E $iTn.(.,.)f/, 

|/3| = |a|-2 

/3' de ^^-^^^ 

d=l /3eN° ■ ij=l,---,|a| 73e {13) 

\l3\ = \a\~2 if^j 

— Z^ 2^ Q,] g^idg^jd g^Ja-e^^-e^,.- 

d=l ij = l,--- ,|a| 

Since G is an orthogonal matrix, one has 

D 
E 9{)^dg-djd = h,-dj ■ (4.52) 

d=l 

Thus ()4.5ip can be further simplified as 

'' = . . 1^ , , «..(«.. -1)5^^"-^-^ = 1^ 9^^"-^-' ^'-''^ 

ij=l,---,|o| »^ » d=l 

which completes the proof of ()4.43p . The equations (|4.44p(|4.45p(|4.46p can be proved using 
exactly the same technique. The detailed proofs are omitted here to avoid redundancy. D 

It is not difficult to find that (j4.42p and (j4.43p still hold if we replace t with Xj or Xj 
for any j = 1, • • • , D. Such observation leads to the following two lemmas: 

Lemma 8. The following equalities hold for arbitrary a € N .■ 

D D r.~ D D . ., r. 

mi-^.^^f^-^^ = Y.Y. E ^^^M(^)^V)^.'^3fo{^)~e,, (4.54) 

j=l d=l ^ j=l d=l ip(zx>W\ ' ^ 

E E 5^^i^"-2ed = E E E ^;^n3(^(a),^) — Uj/,(^)_2e,. (4.55) 

j=i d=i ^ j=i d=i (peDkl ■ ^ 
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Proof. Replacing t with Xj in ()4.42p . we obtain 

D D , ^, ^ D D r.~ 

j=l d=l (^gDl^l ' -^ J=l d=l ^ 

D D D r.~ D D r^~ 

j=l d=l 1=1 * d=l 1=1 * 



(4.56) 



This equation is the same as (j4.54p . The proof of (j4.55p is almost the same. D 

Lemma 9. The following equalities hold for arbitrary a G N^ : 

D D rs- D D , .^ r. 

Yl 2^("i + l)^/"-<=rf+<=. =Y.Y Yl ^^"'?(^(")' '^^ ■ ('^j('^) + ^)-^LM-e,+ej, 

j=i d=i ^ j=i d=i y,ex>|Q| ■ ^ 

EE(«^ + 1)^/-2^.+^.=I;E E ^n,(,(a),^).(a,(^) + l)^/.(^).2e,+e,. 

j=l d=l ^ j=l d=l <^eD|a| ^ 

Proof. Replacing t with Xj and substituting a + ej for a in (|4.42p , one has 

D D 

1 
D D 






V-^V-^, ^, V-^ 0"(y')! TT / / N xC^^id r 

2^2^(a, +1) 2^ — ——-n3(?(a + ej ),(/;) — /,(^)„e, 
X^Srsr S^ (o-(v?) + ei)! .5wd 



(4.57) 



j=l a!=l i=l ip£T>\°'\ 

D D 



J 



= EE E ^n,(.(a),^).(a.(,,) + l)|^/.(,),_,. 

This proves the first equality. The second equality can be similarly proved, and the details 
are omitted. D 

Now the proof of Theorem [3] is given as follows: 

Proof of Theorem\^ Since (ni,--- ,n£)) is a unit vector, we let G = {gij)DxD be an 
orthogonal matrix with its first row as (ni,--- ,n£)). Now we use this matrix as the 
rotation matrix and define w as (j4.36p and (j4.32p . It is obvious that the relation between 
w and w is linear. Therefore, there exists a constant matrix R (see (|4.36p ) depending on 
G such that 

w = Kw, (4.58) 

and R is invertible since w can be obtained from w by applying the rotation matrix G~^. 
Lemma [7H9] have clearly shown that the "rotated equations" 

T(*)^ + E T(*)M,(ti;)||^ = (4.59) 
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can be deduced from ()4.6p by linear operations. Thus there exists a square matrix 'ii{w) 
such that 

U{w)Tiw)^ + ^Il{w)Tiw)M,{w)^ = (4.60) 

j=i -^ 

is identical to ()4.59p . Matching the terms with time derivatives, one finds ii{w) = 
T(w)RT(w)-i. Thus ()i:60D becomes 

T(to)^ + fl T{w)-RM,{w)^ = 0. (4.61) 

Using ()4.40p . the above equation can be rewritten as 

T(*)I^^ + EE5.,T(^)RM,(^)^ = 0. (4.62) 

j=i d=i 

Compared with (|4.59p . one concludes 

D 

J^gijT^ i'w)'RM j{w) = T(-u;)Mi(Rw)R. (4.63) 

i=i 

Multiplying both sides by R"^T('Uj)~^, (j4.7p is attained. Recalling Mi = Am and that 
the first component of the macroscopic velocity after the rotation is li • n (see ()4.32p ). 
the diagonalizability and the eigenvalues of the matrix (|4.7p are naturally obtained using 
Theorem [21 D 

5 Riemann Problem 

Though the regularized moment system ()4.3p is given by the moment expansion up to an 
arbitrary order M thus extremely complex, we can clarify appreciably the structures of the 
elementary waves of this system with Riemann initial value, including the rarefaction wave, 
contact discontinuity and shock wave. Definitely, the structure of the elementary wave is 
fundamental for further investigation into the behavior of the solution of the system. Fur- 
thermore, the solution structure of the Riemann problem is instructional for studying the 
approximate Riemann solver, which is the basis of the numerical methods using Godunov 
type schemes. The analysis below shows that the structure of the elementary wave of the 
Riemann problem is quite natural an extension of that of Euler equations, which indicates 
that the regularized moment system (|4.3p is actually a very reasonable high order moment 
approximation of Boltzmann equation. Following [T5] where the multi-dimensional Euler 
equations are studied, we consider the xi-split, £)-dimensional Riemann problem as below: 

f dw ; dw 

-kt + Amtt- = 0, 
at oxi 

( WL if xi <0, (5-1) 

w{xi,t = 0) = < 

[ wji if xi > 0. 

The Riemann problem with ID velocity space has been studied in [1] in detail. Here we 
focus on the case of L* > 2. 
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Let us first recall the definition of the notations Ajvf, Ajv/, B, B, w, Va^i, Ta.i, v^ and 
Aj^fc in Section [31 In particular, we need the expressions of BB^^ and d"', which read 



BB 



I 








B21 


I 





* 


* 


I 



and for any |a| = M, fc = ai + 1, j = A/z)-i(a), 



(i338T l 



(1339? ) 



where Ij is the j'-th column of the N^ x N^ identity matrix. ra,i, which depends on v^^' 
and Xi^k = ui + Ci^kVO, is the eigenvector of Am for the eigenvalue Aj^^ = ui + Ci^kVO, 
where j = A/z)_i(a), k = M + 1 — \a\. As the first conclusion on the Riemann problem 
()5.ip . we have the following theorem: 



Theorem 4. Each characteristic field of ()5.ip is either genuinely nonlinear or linearly 
degenerate. And one characteristic field is genuinely nonlinear if and only ifv (determined 
by the right eigenvector through (|3.23p ) and the eigenvalue A = ui + CVO satisfy one of 
the following two conditions: 

1. V = v^^\ and C subject to HeM+i{C) = and C 7^ 0; 

2. v = v^^\ j = MD~i{2ek), k G V\{1}, and C subject to HeM~i{C) = and C 7^ 0. 

Proof. Let r denote an eigenvector of A.m with the eigenvalue \ = ui + Cy6 and v is the 
corresponding vector determined by (j3.23p . Since 



^D 



V Dp 



(5.2) 



depends only on p, ui and p2ej/2, d £ V, we have 



V,,,A • r 



• prp + 1 • CV6'rp + 



c^e 



D 



2p 

Vec 



DpVO 



P^P + Yj 



d=2 



^•^'^--/^ 



q2\ ^2 



(5.3) 



liv = v^'^\ then ([339?]l shows v^^'^ = BB"i/i and ireA/+i(C) = 0. From (l338T]l . we 
get 



i;i = 1 and Vj^^^_^(2e2) = ■■■ = Vf^^_^^2eD) 
Thus ()5.3p can be written as 

{D + 1)V9C 



2D 



1 



V,,,A • r 



4L>2 



(C^ + D) . 



(5.4) 



(5.5) 



Hence, 



V^^A • f = 0, if C = 0, 
VioA • f ^ 0, otherwise. 



31 



If V = v^^\ j = Mn^i{2ek) for any k G V\{1}, p39T]l shows v = BB~^ Ij and 



HCM^ 


^i(C) 


= U. l^tom (|;-i.58f|), 


we can 


get 








Tjj = 1 and vi = vi 


= 0, 1-- 


= Md^ 


-i(2ed) 


Then 


(1L3D 


is simphfied as 


V^A 


■ ra,i -- 


^^ 



(5.6) 

Again, we have 

j V^X • f = 0, if C = 0, 

[ Vu,A • f ^ 0, otherwise. 

• Otherwise, (|3.58[p indicates vi = t'A/'i5_i(2efc) = for each k G 'D\{1}. Hence Vi^A-f = 
always holds. 

This completes the proof. D 

This theorem reveals that for each characteristic field, the eigenvalue is constant or 
varies monotonically along the integral curve, resulting in simple wave structures. Below, 
some elementary waves including the rarefaction waves, contact discontinuities and shock 
waves are studied in detail, and the basic relations across these waves are established. 

The analysis below is based on the fact that an eigenvector r of A^ for the eigenvalue 
X = ui + CyO depends only on v and C. With Theorem [5] and the forms of B and BB in 
Lemma m we can divide characteristic fields into three cases: 

Case 1: V = v^^\ and C subject to HeM+iiG) = 0, and C 7^ 0. 

Case 2: V = v^^\ j = 7Vi)-i(2efc) for any k G T>\{1}, and C subject to HcM-iiC) = 
0, and C / 0. 

Case 3: otherwise. 

For convenience, let characteristic field a denote the characteristic field corresponding 
to the eigenvector f^,? for the eigenvalue Xi^k = ui + Ci^kVO with i = ai, A: = M + 1 — |a|. 
Below, the rarefaction waves, contact discontinuities and shock waves will be studied 
respectively. 

5.1 Rarefaction waves 

For the regularized moment system, if two states w and w are connected by a rarefaction 
wave in a genuinely nonlinear field a, then the following two conditions must be met: 

• constancy of the generalised Riemann invariants across the wave, saying the integral 
curve w{C) = (w^i(C)) w;2(C), • • • ,wn{C)) ^^ the iV-dimensional phase space satisfies 



w\C) = r&,^{w), (5.7) 



with i = ai. 

divergence of characteristics 



\,k{w^) < Xi,k{w^). (5.8) 
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, Wj^) in the phase space, the 

D 



Fortunately, for a given point w^ = {p^ ,u\ 

integral curve across w^ can be given. Since P = j^ /.P'ied^ ^^ ^^* P^ = j^ Z^P'^Cd' "^^^ 

d=l d=l 

results are rather tedious, and here the integral curves are only partially given in three 
cases as below: 



If V = v^^>, we have 



c,„kVe. 



' ' «d 



Let r 



2D -I 



p2 



^ ^P2e,/2 



2D 



pO, d £ V\{1}. 



and then we have 



p(C) = /9°exp(C), 



MO = 4 + 



r-1 



exp 



UdiO = ul d = 2,---,D, 



P2e, (C) = P2ei + 



^l - D 



DT 

p(C)=p°exp(rC). 



p'^[exp(rC)-l], 



(5.9a) 

(5.9b) 

(5.9c) 

(5.9d) 
(5.9e) 



Iiv = v'^^\ j = 7VD_i(2efc), k G V\{1}, we have 
rp = 0, , rua =0, deV, rp^,j2 = P^ 
Hence, the integral curve satisfies 

piC) = P\ 

Ud{0 = '^% (i=l, 

P2ei (C) = P2ei : 

/'2C 
P{C) = / exp ( — 



' ^P2<^d 



/2 = o, dev\{k}. 



,D, 



(5.10a) 
(5.10b) 
(5.10c) 

(5.10d) 



Otherwise {v = v^^\ j ^ Mo-ii'^e.k) for any k G V), 



Hence, we have 



rp = ru^ = rp^^j2 =0, d^V. 



m = P', 

P2ei (C) = pI 



:ei' 



ui{0 = u'„ 

piC)=P°- 



(5.11a) 
(5.11b) 
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One can check that ()5.9p . ()5.1ip and ()5.10p satisfy ()5.7p . And an eigenvalue, which satisfies 
dMSD, of Aa/(w(C)) is as 



( 



(t«(C)) = ui(C) + c.,fcVp(C)/p(C) 



exp ( ^C ) - 1 



for i;(i), 



^^? + Q,fc\/^exp T-^") , for ^(J), J = AAz3_i(2efc), A; G P\{1}, 

li? + Cj^fcV^, otherwise. 

It is convenient to verify that Sj_fc(wj(C)) ^ Si^fc(i(;'^) if and only if Cj^fcC ^ Oi and t? and 
Cj^fc satisfy case 1 or case 2. Therefore, if the left state w and the right state w are 
connected by a rarefaction wave and let 'uf' = w^, (|5.8p indicates Si^k{w^) < Sj^fc(itJ^), 
hence Ci^kC > and ^,0^^^ satisfies case 1 or case 2. Therefore, we have that 



for case 1: 
and 



if Q,fe > 0, then uf < nf , p^ < p^, (5.12a) 

if Q,fc < 0, then uf < nf , p^ > p^. (5.12b) 



for case 2: 
and 



ni = u^, d-- 


= 2,--- 


,^, 


if Ci,fc > 0, then u{ 


= nf, 


p^ <p^ 


if Cj,fe < 0, then u[ 


= nf, 


/>p^ 



(5.13a) 
(5.13b) 



5.2 Contact discontinuities 



For a contact discontinuities, ()5.7p is still valid, and the divergence of characteristics is 
replaced by 

^kiw"-) = K^kiw"^). (5.14) 

According to Theorem H] and analysis in Section 15. 1|. the contact discontinuities can be 
founded if and only if v and Cj^^ satisfy case 3. 

• For v^^\ ()5.14p means Cj^fc = 0. Substituting it into ()5.9p . we can get Ud, d G P are 
invariant, while p, p2ei are not (otherwise, (j5.9ep gives us ^ = 0, thus w^ = w^). 

• For v^i\ j = 7Vz)_i(2efc), A; G P\{1}, (1^1^ means Q^fc = again. (fUTH shows p, 
Ud, d G T), p2ei are invariant, while p is not (otherwise, ()5.10dp gives us C = 0, which 
results w^ = w^). 

• Otherwise, (|5.1ip shows ui, p, p2ei are all invariant. 

Summarizing the discussion above, we conclude that if Cj^fc ^ 0, then ui, p, p2ei are 
invariant across the contact discontinuities, while if Cj^^ = 0, ui is invariant and p is not. 
However, Ud, d = 2,- ■ ■ ,D may change discontinuously across a contact discontinuity. In 
fact, the case v = B"^/^, d = 2,- ■ ■ ,D corresponds to a contact discontinuity where Ud is 
discontinuous. This is similar as the Euler equations. 
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5.3 Shock waves 

The discussion of the shock wave needs some more scrupulosity. As is well known, the 
jump condition on the shock wave is sensitive to the form of the hyperbolic equations. 
Thus, before we give the Rankine-Hugoniot condition, it is necessary to rewrite (|5.ip in 
an appropriate form. However, (|5.ip cannot be written as conservation laws due to the 
presence of TZ\j jj{a). Nevertheless, (j5.ip can still keep the conservation of the conservative 
moments with orders from to M — 1. Therefore, ()5.ip can be reformulated by MdUM — 
l)e/)) conservation laws and N — MdUM — l)e[)) non-conservative equations. 
Let 

F=iFo,Fe„Fe„--- ,FMen), ^« = A / C f ^t \a\ < M, (5.15) 

where ^" = 0^=1 ^d"^ ^^^'^ -^o stands for Fa\a=o- Then ()5.ip can be written as 



dt dxi 

The relation between F and w is 



^— + (ai + l)^ = 0, a < M, 

ot ox I 

" + (ai + l)^-7ei,^(a) = 0, \a\=M. 



(5.16) 



He, ■ " 



y- [_l)W-P\p^ ! V V^/ g|a^g|/2 



\P\<W\ 



(«-/5)! ' (5.17) 



F F^ 



where Hca I —/= ) = 11^=1 ^^a^ I ~~/= ) ^iid Hea{x) = if at least one aj is negative. In 
addition, F^ is as 

Hea+ei-l3 I -7= 

For convenience, the quasi-linear form of ()5.16p is written as 

dF dF 

where r(i^) is an A^ x A^ matrix and depends on ()5.16p . 

Since (j5.19p is not a conservative system, we have to adopt the DLM theory [lUj to 
study the shock wave. For a shock wave the two constant states F and F are connected 
through a single jump discontinuity in a genuinely non-linear field a travelling at the speed 
Sa, and the following two conditions apply 

• Generalized Rankine-Hugoniot condition: 

J [SaI-T{^{i.;F^,F^))]^{u;F\F^)du = 0, (5.20) 
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where I is the NxN identity matrix, and ^{v^F , F ) is a locahy Lipschitz mapping 
satisfying 

*(0; F^, F^) = F^, *(1; F^, F^) = F^. (5.21) 

We refer the readers to [lOj for details. 

• Entropy condition: 

Ai,fc(F^) >Sa> Ai,fc(F^), (5.22) 

where i = ai and k = M + 1 — \a\. 

For conservation laws, ()5.20p is the same as the classical Rankine-Hugoniot condition. 
Thus the first MdUM — l)e£)) rows of (|5.20p are independent of $. This allows us to 
analyze the properties of the shock waves without regarding the form of $. 
The first equation and the {D + l)-th equation of ()5.20p are precisely as 

Af-pV = ^a(p^-p''), (5.23) 

P''inif+Pk-P''ini?-P2e, = 5«(p^nf - Af). (5.24) 

• If p^ / p^, ([5:23]) and (fOil) give 

o u o u 

Dq ^ f p — (5.25a) 

pL _ pH 



p^{u{f +pg^ - p^{uff-pg^ 



p^u{ — p^uf 



(5.25b) 



Substituting ()5.25ap into (j5.22p and multiplying both sides with {p^ — p^)"^, we get 

/(nf - nf )(p^ - p^) > Q,fc(/ - pyVeR, (5.26a) 

p'^iuf - ^f )(/ - p^) < Q,fc(/ - p^Ve^. (5.26b) 

If Ci fc > 0, dOeil) gives 

(uf-nf)(/-p«)>0. (5.27) 

Therefore, we can divide (|5.26p by {u^ — u^){p^ — p^) and obtain 

/ , Q,,(p^-P^) ^ p^ 

Thus we have 

{p^fe^ > {p^fe^. (5.29) 

Furthermore, (|5.25p gives us the relation 

(P^^ - p''){pL, -pD = P''p\u\ - nf )^ (5.30) 

If p^ < p^, (j5.27p indicates uf < uf . (|5.29p can be written as p^p^ > p^p^, so we 
have p^ > p^. If p^ > p^, (|5.27p indicates uf > u^. Summarizing these results, we 
get 

if p^ < p^, then nf < uf and p^ > p^, 

r R r R "^ "^ ' (5.31) 

if p^ > p^, then uf > uf. 
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Analogously, if Cj^fc < 0, we have 

K--f)(/-/)<o, (pY^^<(/)^0^ 

and ()5.30p still holds. Hence we get that 

if p^ > p^, then uf < uf and p^ <p^, 
if p^ < p^, then nf > uf . 



(5.32) 



(5.33) 



If p^ = p^, (j5.23p and ()5.24p make that uf = u^ and P2e ~ P2e ' respectively. 
Therefore, (|5.22p is turned into 



The following result is then attained 



R ^L ^ R 



if Ci,fc > 0, then Ui = nf , p^^ > p". 



fl „L ^ R 



if Cj^fc < 0, then uf = nf , p^ < p 



(5.34) 



(5.35a) 
(5.35b) 



Now we summarize all our discussion on the entropy conditions of the three types of 
elementary waves in the following theorem: 

Theorem 5. For the Riemann problem ()5.ip . for the wave of the a-th family, Ci^k, the 
macroscopic velocities and pressures on both sides of the wave have the relation with the 
type of the wave as in Tablel^ where Ci^k corresponds to the eigenvalue Xi^k = "Ui + Cj.fcVP, 
i = ai and A: = M + 1 — |a|. 



Wave type 


Eigenvalue 


Velocity and Pressure 


Rarefaction wave 


Q,fc > 


u{ < lif , p^ < p^ 


Ci,fc < 


u{ <u^,p^> p^ 


Shock wave 


Q,fc > 


uf < nf , p^ > p^ 


Q,fc < 


u[ < nf , p^ < p^ 


Q,A: 7^ 


u[ > nf 


Contact discontinuity 


Q,fc = 


u{ = nf 


Ci,fc ^ 


nf = nf , p^ = p^ 



Table 1: The relation between the type classification of elementary wave and the eigen- 
value, macroscopic velocity and pressure. 
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Appendix 

A Collection of Notations 

We list below some of the notations used for convenience. 

Q{f, f) The Boltzmann collision operator taken as zero in this paper. 

Ci^k The i-th zero of Hek{x). 

a D-dimensional multi-index. 

a D-dimensional multi-index, defined as d = a — aici. 

a {D — l)-dimensional multi-index, defined by ()3.6bp . 

Sd,m All the D-dimensional multi-index satisfying the sum of components 

not more than M, defined by ()3.6cp . 

SD,M{<y) The subset of Sd,m, defined by ()3.6cp . 

Nd{oi) The ordinal number of a in Sd,m, which is permuted by lexicographic 

order, defined by ()3.7p . 

N Cardinality of Sd,m- 

w The basic moments in the moment system, defined by ()3.9p . 

Am The coefficient matrix of the moment system, defined by (|3.10p . 

I Identity matrix. 

TZm,d Regularization term based on the characteristic speed correction. 

TZ\,j jy Regularization term in Xj direction. 

Am The regularized matrix of Am, defined by ()3.19p . 

A The diagonal matrix, used to make Am dimensionless, defined by 

Am The dimensionless matrix of Am — uil. 

Pei+sk ) 9a The dimensionless of Pa+Ck i fa- 

r The right eigenvector of Aj\,./ , defined by p.30p . 

V The components of it is a subset of that of r, and vj\fjy_^(^^-j = r^j^ .g^-. ■ 

B The coefficient matrix defined in Lemma El 

B The block diagonal matrix, used to make B a unit lower triangular 

matrix, defined in Lemma [2j 

«(■'■) The basis of v, defined by (13391) . 

r^^j The eigenvector of Am, depended on v^^\ j = A/'D-i(a), for the eigen- 

value Cj^fc, k = M + 1 — \a\. 

Vd,m The characteristic polynomial of Am, defined by (|3.63p . 

\,k The eigenvalue of Am, h,k = ui + Ci±\/9. 

r^^i The eigenvector of Am, corresponding to r^, j, for the eigenvalue Cj^fc, 

k = M + l-\a\, defined by (IMjl . 

Vd,m The characteristic polynomial of Am, defined by ()3.67p . 
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G = {gij)DxD The rotational matrix. 

w The rotated moments, defined by ()4.32|) and (|4.36p . 

iD, ip The Grad-type indices in V^. 

ij(-) Conversion from multi-indices to Grad-type indices. 

o"(-) Conversion from Grad-type indices to multi-indices. 

^g{'&-, (p) Product of entries in the rotational matrix, defined by (j4.24|) . 

Al^ The set of all m-permutations in {1, • • • , n}, defined by ()4.2ip . 

T{a) The set of all Grad-type indices whose corresponding multi- indices are 

a, defined by (fil^ . 
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